arXiv:1501.00459vl [hep-ph] 2 Jan 2015 


EPJ Web of Conferences will be set by the publisher 
DOI: will be set by the publisher 

© Owned by the authors, published by EDP Sciences, 2015 


Nucleosynthesis constraints on the faint vector portal 

Anthony Fradette 1 , Maxim Pospelov 1,2 , Josef Pradler 3 3 , and Adam Ritz 1 

1 Department of Physics and Astronomy University of Victoria, Victoria, BC V8P 5C2, Canada 

2 Perimeter institute for Theoretical Physics, Waterloo, ON N2J 2W9, Canada 

3 Institute of High Energy Physics, Austrian Academy of Sciences, A-1050 Vienna, Austria 


Abstract. New Abelian (7(1)' gauge bosons V t , can couple to the Standard Model through mixing of the 
associated field strength tensor V jlv with the one from hypercharge, F^ v . Here we consider early Universe 
sensitivity to this vector portal and show that the effective mixing parameter with the photon, k, is being probed 
for vector masses in the GeV ballpark down to values 1CT 10 < k < 1(T 14 where no terrestrial probes exist. The 
ensuing constraints are based on a detailed calculation of the vector relic abundance and an in-depth analysis of 
relevant nucleosynthesis processes. 


1 Introduction 

The origins of our Universe may well be rooted in infla¬ 
tion or alternative cataclysmic scenarios that regard the 
very earliest moments of existence. However, despite the 
impressive success of observational cosmology over the 
past decades, the earliest true direct window into the be¬ 
ginnings remain observations of light element abundances. 
They concern the epoch of primoridal nuclear transforma¬ 
tions at cosmic times r > Is. The overall concordance 
of the Big Bang nucleosynthesis (BBN) predictions with 
the observationally inferred primordial values is one of the 
most impressive successes of modern day cosmology and 
particle physics. Today, BBN is used as a toolbox to put 
models of new physics to a stringent test Jlj], whenever 
they predict some interference with the the standard pro¬ 
cesses in the observable sector at t > Is. 

Under the assumption of a canonical sequence of cos¬ 
mological events, the Universe emerged from inflation and 
baryogenesis much prior to BBN. Such sequence then al¬ 
lows one to put stringent constraints on very weakly inter¬ 
acting sectors of new physics beyond the Standard Model 
(SM). The kinetic mixing of a new U{\)' vector V M with 
hypercharge Fj v V /iv is of particular interest as the mix¬ 
ing with the photon leads to numerous experimental con¬ 
sequences and much attention was devoted to this vector 
portal in recent years j2l. Below the electroweak scale, the 
coupling of V to the SM is essentially given by its mixing 
with the photon 0], 

£v = - K -F^ = eKVXv (1) 

With k and my being the only free parameters, the model 
provides a simple, and technically natural prototype sce¬ 
nario for a light, weakly interacting new particle sector. In 
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the following we will concentrate on a Stiickelberg ori¬ 
gin of my that allows to maintain gauge invariance in 
U{\)' without complicating the phenomenology by hidden 
Higgs particles IT] see e.g. |4] for the phenomenology and 
|f5ll for cosmological constraints on the latter scenario. 

The SM decay modes of V are well known. When 
hadronic decays are kinematically accessible, one can use 
experimental data on the A’-ratio to infer couplings to pho¬ 
tons in the time-like direction, and hence to determine the 
decay rate Ty and all branching ratios. Below the di-muon 
threshold and for my > 1 MeV the vector V decays to 
electron-positron pairs only, thereby setting its principal 
lifetime. 


Ty 


3 

k 2 aniy 


= 270 s x 


1 GeV / 10~ 12 \ 2 

my \ K ) 


( 2 ) 


where a is electromagnetic fine structure constant. In the 
following the cosmological consequences of U{\)' vectors 
with masses in the MeV-GeV range, and lifetimes long 
enough for the decay products to directly influence primor¬ 
dial nucleosynthesis are explored. These vectors have a 
parametrically small coupling to the electromagnetic cur¬ 
rent, and thus an extremely small production cross section 
for e + e~ —> Vy, cr pro d ~ K 2 na 2 /s ~ 10~ 54 cm 2 where we 
took = 200 MeV and k = 10~ 12 from above. Such 
small couplings render these vector states completely un¬ 
detectable in terrestrial particle physics experiments. 

Despite the tiny production cross section, any charged 
SM state that is populated in appreciable number in the 
early Universe at temperature T ~ my may yet emit V. 
With the above ballpark numbers in (0, parametric es¬ 
timates suggest that an amount of MeV/baryon may be 
stored in V-particles. Followed by late decays back to SM 
states, visible energy is therefore being injected into the 
primordial plasma at levels that are probed by BBN. The 
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Figure 1. Dark photon inverse decays are the leading contribu¬ 
tion to sub-Hubble production rates in the calculation of the V 
relic density. 


early Universe is therefore likely to be the only “labora¬ 
tory” where such very dark photons with k ~ 1 0” 12 and 
smaller are being probed. Here we report on the detailed 
analysis performed in [0 where also CMB limits on even 
later decays were considered. Previous partial discussions 
of cosmological signatures of decaying dark photons may 
also be found in 00. 

In the following we assume no other light states x that 
are charged under U(l)'. Therefore, there are no decays 
V —> XX that potentially drain visible modes and thereby 
ameliorate the derived limits from BBN. Using some re¬ 
cent insight about the in-medium production of dark vec¬ 
tors fid ( see a l so 110]) we first discuss the production of 
dark vectors in the next section and explore constraints on 
V-decays into SM in Sec.[3]before concluding with Sec. [4] 

2 Abundance prior to decay 

The relic abundance of weakly coupled dark photons prior 
to their decay is obtained through a calculation of the leak¬ 
age from the observable sector to the hidden sector with 
sub-Hubble rates, r pro d/// <sc 1. This “freeze-in” process 
is dominated by inverse decays of V , through the coales¬ 
cence of e ± , fJ ± ..., and, similarly, through hadronic con¬ 
tributions; see the illustration in Fig. [T] 

The Boltzmann equation for the total number density 
of V takes the form 


abundance can be performed explicitly, 

9 m\T v ^e 
v An ( Hs) T=mv ' 

While further leptonic production channels are easy 
to be included into (0), hadronic production channels re¬ 
quire assumptions about the primordial hadron gas and 
the strength of interaction with photons, denoted by ar e ff. 
At temperatures above the QCD confinement scale T c ~ 
200 MeV light quarks are deconfined and individual quark 
contributions can be added to Y v in a straightforward man¬ 
ner. Below T c one may use a free gas of mesons as an ap¬ 
proximation to the hadronic (non-baryonic) particle con¬ 
tent in the early Universe. The production via inverse 
charged pion and kaon decays \n'n ,KK) —> V can 
then be included using a scalar QED model with effective 
coupling strengths like o^g(my) — K 2 a* n ( sfs = my) where 
a nn is extracted from BaBar cross section measurements 
of e + e~ —> y* —* n + n~(y) JH , and similarly for charged 
kaons fl0 . 

Finally, there is a possibility of resonant production of 
V by virtue of the thermal bath. Such in-medium effects 
may be cast into an effective mixing angle. 


«t.l 


k 

|i-n r ^/<r 


( 6 ) 


with 1 \jj being the transverse (T) and longitudinal (L) 
photon polarization functions in the primordial, isotropic 
plasma. The expressions for Y\ jj, can e.g. be found in 11311 ; 
the longitudinal polarization function [0 used here is, 
n£ ere = ml Her - m 2 )nf f ^ and e> is the (dark) pho¬ 
ton energy. Equation (0 informs us about the condition of 
resonant dark photon production. 


Re n r.ziw, T r _ t ,l) = my- (7) 


The condition depends on temperature T as FU / are pro¬ 
portional to the plasma frequency, ejp(T). Most impor¬ 
tantly, the resonance temperature T r jjXoj) as a function of 
frequency e) is parametrically larger than my with a mini¬ 
mum frequency at which the resonance can happen. 


hy + 3 Hny = C. (3) 

The right hand side is the collision integral and provided 
that V never reaches thermal equilibrium, it gives the num¬ 
ber of V states emitted per unit volume and unit time. In 
the Maxwell-Boltzmann approximation for the SM distri¬ 
bution functions, and in the limit that only electrons coa¬ 
lesce, the integration can be done analytically 001 (see 
also 0) 

C * ^Yy^mlTKHmy/T). (4) 

Here, r= K 2 amy/ 3 is the decay width of V to elec¬ 
trons, up to corrections ( m e /my) 2 and K\ is the modified 
Bessel function of the second kind. In terms of the number 
of V particles normalized to entropy density, Y v - ny/s, 
the cosmic time integral in (0 for the final “freeze-in” 


7|\rnin — my 


2na 


1/2 


^ 8 my. 


( 8 ) 


Thus resonances occur at parametrically larger tempera¬ 
tures (by a -1 / 2 ) than my, for which H(T) is significantly 
larger than at T =* my at which the V freeze-in production 
has its biggest contribution. Therefore, resonant contri¬ 
butions to Yy do not alter the picture drastically though 
numerically they may constitute as much as 30%. 

After production, the momentum of V redhifts quickly 
so that at the time of decay the energy of V is to good ap¬ 
proximation given by the rest mass, Ey = my. The decay 
deposits this energy into leptons, hadrons, and hadronic 
resonances. The energy prior to decay that is stored per 
baryon is therefore given by 

£p.b. = myYy —, (9) 

fib. o 
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where «/,,o/ s o = 0.9 x 10~ 10 is the baryon-to-entropy ra¬ 
tio today. Equipped with £p.b. as a function of m v and 
k following the detailed calculations of the V “freeze in” 
abundance in 0], we may now explore its consequences 
for BBN. 

3 F-decays during BBN 

Primordial nucleosynthesis predictions are affected for 
dark photon decays with cosmic lifetime t > 1 s or larger. 
Ensuing constraints are then governed by a combination of 
lifetime and abundance, both being complementary with 
respect to the vector mass: Ty ( Yy ) decreases (increases) 
with larger my. From this one expects constraints as local¬ 
ized islands in those parameters where the epoch of pri¬ 
mordial nucleosynthesis exhibits its greatest sensitivity. 

3.1 Major effects and treatment 

The effects on BBN are understood by considering electro¬ 
magnetic and hadronic energy injection separately. Prior 
to decay, the V abundance relative to baryons is substan¬ 
tial, ny/nt, < 10 8 for Ty < 1 s, and the decays of V in¬ 
ject electrons, muons, and mesons in numbers larger than 
baryons. 

Dark photon decays with my < 2m„± = 279 MeV re¬ 
sult exclusively in injection of electromagnetic energy, be¬ 
cause V —» e + e~,p + p~ are the only kinematically accessi¬ 
ble modes. Muons typically decay before interacting, and 
electron-positron pairs are quickly thermalized by interac¬ 
tions with background photons. The resulting electromag¬ 
netic cascade with spectrum f y (E y ) entails a large number 
of non-thermal photons that may then spall light elements. 

Importantly, the spectrum has a relatively sharp cut-off 
for energies above the e ± pair-creation threshold, E p;ur - 
ml/(22T). Photons with E y > Ep a j r are being dissipated 
before they interact with nuclei, and to good approxima¬ 
tion f y (E y ) = 0 for E y > £p a i r . Photons with E y < E pa i r , 
however, undergo slower degradation processes and may 
interact with the light elements before being thermalized. 
Equating E pa ; r against the photo-destruction thresholds (in 
brackets below) yields the temperature and thereby the 
cosmic time f p h of biggest impact for a spallation channel: 

( 2 x 10 4 s, 7 Be + y -> 3 He + 4 He (1.59 MeV), 

fph^i 5 X 10 4 s, D + y->n + p (2.22 MeV), 

( 4 x 10 6 s, 4 He + y —> 3 He/T + n/p (20 MeV), 

The spallation rate of species N with number density 
lift is given by 

i~*E max 

r ph (r) = 2«jv| dE y f y (E y )o- y+N ^x(Ey), (10) 
d /',h- 

where <r y+ x-,x(Ey) is the photo-dissociation cross section 
for y + N —> X with threshold £ thr . The factor of two 
accounts for the two independent cascades that form in a 
back-to-back decay of V at rest, each with a maximum 
energy of £ max = max|E pa i r , E; n j/2j. All spallation reac¬ 
tions listed in Cl are taken into account in the numer¬ 
ical analysis. We note in passing that neutrino injection 



Figure 2. Pythia simulation on the average number of particles 
per V decay and on the injected electromagnetic energy after all 
particles have decayed or annihilated to electrons and photons. 
About one third of the energy is carried away by neutrinos. Nar¬ 
row hadronic resonances are neglected. 


from muon decay constitute only minor corrections to the 
photon-induced processes listed above d. 

For vector masses above the di-pion threshold, my > 
2m n ±, hadronic modes are accessible in the decay of V and 
the effects on BBN are more intricate. In the hadronic de¬ 
cay of V only n ± , K ± , and /f/ , with lifetimes r ~ 10 8 s, 
and (anti-)nucleons have a chance to undergo a strong in¬ 
teraction reaction before decaying by themselves. 

Before deuterium formation at T ^ 100 keV, only 
charge exchange reactions on nucleons, such as n~ + p — » 
j r° + n, are possible. They change the n/p ratio and thereby 
most prominently the primordial 4 He value. After the deu¬ 
terium bottleneck—once light elements have formed— 
charge exchange creates “extra neutrons” on top of the 
residual and declining neutron abundance. In addition, 
absorption with subsequent destruction of light elements 
such as 7r~ + 4 He —> T + n is now operative. Spallation of 
4 He may also have a secondary consequence: the produc¬ 
tion of mass-3 elements with non-thermal kinetic energy 
may induce reactions of the sort T + 4 Heb g —> 6 Li + n. 
In the numerical analysis, these processes as well as sec¬ 
ondary populations of n L from kaon decays, and hyperon 
producing channels from reactions of kaons on nucleons 
and nuclei are being accounted for. Furthermore, in our 
analysis, we restrict ourselves to reactions at threshold, 
with charged pions and kaons being thermalized before re¬ 
acting on light elements; such approximation generally re¬ 
sults in more conservative constraints. A detailed quantita¬ 
tive discussion of incomplete stopping can be found in @1- 

Finally, baryon pairs are produced in the V-decay for 
my > 2GeV. Final state nucleons n and p will prefer¬ 
entially annihilate on protons with an annihilation cross 
section (cr ann u) ~ . The injection of nh then results in 
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Figure 3. Final state branching ratios of long-lived mesons and 
other relevant decay products. BaBar measurements of the e* —» 
n* and e* —> K* cross sections up to my = 1.8 GeV are stiched 
together with a Pythia simulation starting at my > 2.5 GeV; 
a branching to K L was neglected and the fraction of m v that 
is ultimately being converted to electromagnetic energy is la¬ 
beled Br cm . 


one net p —* n conversion with associated energy injection 
of m p + m n . Annihilation on neutrons with similar cross 
section is also possible and pp injection results in one net 
n —> p conversion. Assuming equal cross sections, the 
relative efficiency is p/(n + p) and n/(n + p), respectively. 

At threshold, the rate for neutron injection can be in¬ 
ferred from a measurement of the e + e~ —> nh cross sec¬ 
tion, CT e + e -—*nn ~ 1 nb d. Given a total hadronic cross 
section cr c+e -^ had ~ 50 nb at this energy, the branching 
fraction to a neutron-antineutron pair is ~ 2%. Away from 
the di-nucleon threshold, with multi-pion(kaon) produc¬ 
tion and decays to hyperons and baryonic resonances be¬ 
ing prevalent, V-decays may be simulated using Pythia. 
The ultimate yield of n ± ,K ± , K L , and nucleons prior to 
their decay is shown in Fig. 0 dots depict the average 
electromagnetic energy that is injected after all particles 
have decayed to electrons and photons; e + have been an¬ 
nihilated on e~. The rest of the decay-energy is carried 
away by neutrinos. At lower energies, decay events are 
eventually dominated by two body decays. Above the di- 
pion (di-kaon) threshold, we use BaBar measurements of 
the e ± —> 7 t* and e ± —> K* 1 cross section until an reported 
energy of \[s — my = 1.8 GeV. Relevant ultimate branch¬ 
ing ratios are shown in Fig. 0 the effects of K t are, for 
simplicity, neglected in our BBN anaslysis. 


A more detailed discussion along with a list of all in¬ 
cluded reactions can be found in the original paper as 
well as in the preceding work Q]. Numerical results were 
obtained by usage of a Boltzmann code that is based on 
Ref. 111611 . with significant improvements and updates as 
detailed in ||5t]. Standard BBN yields are found to be in 


agreement with d when using a baryon asymmetry of 
r/b = 6.2 x 10~ 10 and a neutron lifetime of t„ = 885.7 s. 


3.2 Light element observations 

BBN sensitivity is attained by the observational inferrence 
of light element abundances and their estimated error bar. 
Here we briefly discuss those observations that form the 
basis of our obtained regions of interest. 

The most abundant element after hydrogen is helium. 
Its mass fraction Y p is inferred from extragalactic HII re¬ 
gions, and values in the range 


0.24 < Y p < 0.26 


( 11 ) 


have been reported over the years. Owing to potential sys¬ 
tematic uncertainties 018L119i l we adopt (fTTb as the cosmo- 
logically viable range. 

Among recent developments, the precision determina¬ 
tion of D/H from high redshift quasar absorption systems 
stands out itiblEHl. Error bars have reduced by a factor of 
five in comparison to previously available determinations. 
The weighted mean now reads |_2l|]. 


D/H = (2.53 + 0.04) x 10 


-5 


( 12 ) 


D astration on dust grains is, however, a potential source 
of systematic uncertainty, and values as high as 4 x 10~ 5 
have also been reported 1221 23]. In light of this, we adopt 
an upper limit of, 


D/H < 3 x 10~ 5 . (13) 


as well. Finally, producing too little D/H yields a robust 
limit because no known astrophysical sources of this frag¬ 
ile light element exist. We therefore either use the nominal 
lower 2cr-limit from (fl2T> or require (robustly), 

3 He/D <1 (14) 


instead. The latter value is derived form solar system ob¬ 
servations |24]. 

Finally, and with much smaller abundance, the pri¬ 
mordial value of 7 Li/H [25], is lower than the lithium 
yield from standard BBN by a factor of 3-5, 7 Li/H = 
(5.24/y 7 ')x 10“ 10 0. We consider lithium being in con¬ 
cordance with observations if BBN predictions yield 


1(T 10 < 7 Li/H < 2.5 x 1(T 10 . 


(15) 


While new physics may be at the heart of the lithium prob¬ 
lem, we caution that astrophysical depletion mechanisms 
may also play their part in solution to this long-standing 
puzzle; see [26]] for a recent review. 


3.3 Results 

Our results from V-decays and their effect on BBN are 
presented in the my, k parameter space in Fig. [4] Contours 
of constant lifetime, t v and freeze-in abundance ny/ni, are 
shown by the diagonal solid and dotted lines, respectively. 
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Figure 4. Vector mass m v and kinetic mixing parameter k pa¬ 
rameter space with BBN sensitivity. Diagonal gray contours 
depict t v (solid) or n v /n b prior to decay (dotted). Shaded re¬ 
gions are excluded by observations as labeled. The solid (orange) 
closed line is a 2cr constraint from underproduction of D/H de¬ 
rived from G3- Dashed black lines show decreasing levels of 
7 Li/H, 4 x 1CT 10 and 3 x 10 -10 , from outer cirlces to the inner 
ones, respectively. Along the dotted line 6 Li/H = 10~ 12 signify¬ 
ing an two orders magnitude enhanced 6 Li yield, that is, however, 
not yet constrained by observations. 


The regions labeled I-III are in conflict with observations 
as detailed in the previous section. 

In regions I, V decays to e + e~ result in electromagnetic 
energy injection. Region la (t v ~ 10 5 s) is marked by a de¬ 
struction of 7 Be and D. The 7 Li/H abundance is reduced 
to 4 x 10 10 and 3 x 10 10 from the outside to the inside, 
respectively. However, cosmologically favored smaller 
7 Li/H abundances are challenged by 3 He/D < 1 (pink 
shaded region). Using (IT2l ). lower 7 Li/H values are ex¬ 
cluded by the nominal 2 cr lower limit on D/H as depicted 
by the solid closed line. Region lb is additionally marked 
by spallation of 7 Li and 7 Be from non-thermal photons. 
This results in direct production of 6 Li/H > 10 17 —values 
yet too low for being observationally constrained at the 
moment. Finally, Region Ic with Ty ~ 10 7 s is marked 
by 4 He dissociation and net creation of 3 He/D ruling out 
this parameter region. Secondary production of 6 Li is not 
efficient enough to yield an additional limit. 

In region II, Ty < 100 s and V decays before the 
end of the D-bottleneck (T ~ 100 keV). Injection of pi- 
ons, kaons, and nucleons, results in anomalous n p 
inter-conversion. The consequence is an elevated n/p- 
ratio and therefore enhanced D and 4 He yields. The low- 
lifetime/high-abundance region II is correspondingly dis¬ 
favored by Y p < 0.26 and D/H < 3 x KU 5 . 

Finally, region III is marked by the production of “ex¬ 
tra neutrons” at t ~ 10 3 s from V —> nn and from charge 
exchange of n on protons, tt p —> nn° or n~p —> ny. 
In addition, hyperon production by “s-quark” exchange of 
K~ on protons may also result in extra neutrons. With it 
comes a path that may deplete lithium, 7 Be + n —> 7 Li + p, 


followed by 7 Li+/; —> 4 He+ 4 He. With a reduced Coulomb 
barrier, 7 Li is more susceptible to proton burning in the 
second step and the declining 7 Li trend is depicted by the 
dashed lines in Fig. [4] Most of the extra neutrons, how¬ 
ever, end up being captured by protons and the associated 
D/H constraint (fl3l > is given by the orange region. 


4 Conclusions 


The kinetic mixing of a new U(l)' gauge group with the 
Standard Model (7(1) factors of hypercharge and, below 
the electroweak scale, of electromagnetism is one of the 
few portals to the hidden sector with renormalizable cou¬ 
plings. The associated gauge boson V is often called a 
“dark photon” and in this manuscript we have reported the 
cosmological limits from BBN as they have been derived 
in 0]. BBN sensitivity reaches photon kinetic mixing pa¬ 
rameters of k ~ 10 _ 1 4 for 1 MeV < m v < lOGeV, unchal¬ 
lenged from terrestrial dark photon searches, see, e.g. the 
works and presentations [27-33] 

The presented limits are based on a thermal abun¬ 
dance of V and a standard cosmological history of 
the Universe— i.e. uneventful until V-decay—with reheat 
temperatures in excess of my. Additional contributions to 
the V-abundance such as from an initial V-condensate af¬ 
ter inflation may only strengthen the derived bounds. The 
latter source of primordial V-particles is particularly in¬ 
teresting in the context of smaller V-masses. Below the 
di-electron threshold (not considered in this work), V has 
a naturally long lifetime with V —> 3y being the only de¬ 
cay mode. Therefore V can even be a dark matter candi¬ 
date IS0 and ensuing constraints on the photoelectric 
absorption of V on atoms in dark matter detectors have 
started to receive attention only very recently [ 34 - 33] ■ 
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